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Abstract. We analyze quantal Brownian motion in d dimensions using the 
unified model for diffusion localization and dissipation, and Feynman- Vernon 
formalism. At high temperatures the propagator possess a Markovian property 
and we can write down an equivalent Master equation. Unlike the case of the 
Zwanzig-Caldeira-Leggett model, genuine quantum mechanical effects manifest 
themselves due to the disordered nature of the environment. Using Wigner picture 
of the dynamics we distinguish between two different mechanisms for destruction 
of coherence: scattering perturbative mechanism and smearing non-perturbative 
mechanism. The analysis of dephasing is extended to the low temperature regime 
by using a semiclassical strategy. Various results are derived for ballistic, chaotic, 
diffusive, both ergodic and non-ergodic motion. We also analyze loss of coherence 
at the limit of zero temperature and clarify the limitations of the semiclassical 
approach. The condition for having coherent effect due to scattering by low- 
frequency fluctuations is also pointed out. It is interesting that the dephasing 
rate can be either larger or smaller than the dissipation rate, depending on the 
physical circumstances. 



1. Introduction 

Classical Brownian motion is described by the Langevin equation 

mx + ^ (1) 

where x is the position of the particle, rj is the friction coefficient, and T = —WU (x, t) 
is a stochastic force. This equation is meaningful only in a statistical sense. The time 
evolution of a phase space distribution p(x, p) is obtained by solving (|^) for various 
realizations of the stochastic potential, and then averaging over all these realizations. 
The stochastic potential is zero on the average and its correlations are 

{U[^' ,t")U{-K' , t')) = c^{t"-t') ■ w(x"-x') (2) 

Typically these correlations are characterized by a 'short' temporal scale Tc and a 
'microscopic' spatial scale i. Usually, it is further assumed that higher moments 
are determined by Gaussian statistics. The Langevin description can be derived by 
considering a general Hamiltonian of the form 

U = UQ{yi,Tp) +Uenv{'^,Qa-,Pa) (3) 

where x and p are the canonical variables that correspond to the distinguished 
degree of freedom, TLq = Tp^ /2m is the free motion Hamiltonian, and {Qa,Pa) are 
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environmental degrees of freedom. The dynamical variable x may represent the 
position of a large particle. In the one-dimensional version of Brownian motion it 
may represent the position of a piston. The actual conditions for having the reduced 
Langevin description turn out to be quite weak The environment should consist 
of at least 3 degrees of freedom with fast chaotic dynamics. Fast implies that the 
classical motion is characterized by a continuous spectrum with high frequency cutoff, 
such that the motion of the environment can be treated adiabatically with respect to 
the slow motion of the particle. It is essential to assume that the following condition 
is fulfilled 

V 1 

Generic Brownian Motion <^ — ^ — (4) 

£ Tc 

where Tc and £ characterize the correlations of the stochastic potential which is 
experienced by the particle. Eq.(^) is the condition for using the white noise 
approximation (WNA). namely, the noise is characterized by its intensity 

(l>(T)dT ■ |w"(0)| (5) 

J 

and its temporal correlations are adequately described by the formal expression 
0(t) = i/(5(t). Without loss of generality we assume from now on the normalization 
w"{Q) = — 1. The correlations of the stochastic force satisfy 

{mHt'))at X - Ht-t') (6) 

In the general case, (|^) is less informative than (|^). However, in case of a classical 
particle that experiences white noise, the additional information is not required at all! 
It should be emphasized that this observation does not hold if (||) is not satisfied: If 
Tc is larger than £/v, then the particle will perform a stochastic motion that depends 
crucially on the "topography" of the stochastic potential. 

The classical analysis further reveals 0] that associated with the stochastic 
potential there is also a dissipation effect. If the environment is characterized by 
either micro-canonical or canonical temperature T, then the friction parameter will 
be 

77 = — Fluctuations Dissipation (7) 

2kBT 

Thus, any generic environment, in the sense specified above, leads to the universal 
ohmic behavior. The motion of the particle is determined by the interplay between 
the friction and the noise. The friction leads to damping of the particle's velocity, 
while the noise pumps energy back into its motion. Eventually we have diffusion with 
coefficient D^j = i^/v^- If (jj) is not satisfied, then we will have diffusion even in 
the absence of dissipation. The latter " non-dissipative diffusion" is characterized by 
the coefficient Dq ^ £v, where v is the velocity of the particle. This latter type of 
diffusive behavior should not be confused with the generic " dissipative diffusion" that 
characterizes Brownian motion. 

One wonders whether there is a well defined quantized version of the above 
Langevin equation. Is it possible, quantum-mechanically, to characterize a universal 
dynamical behavior that corresponds to classical Brownian motion? Are stronger 
physical conditions required in order to guarantee generic behavior? In the present 
paper we have no intention to give a full answer to all these questions. Rather, as 
in previous publications we follow the Caldeira-Leggett strategy Namely, we 
consider the motion of a particle under the influence of an effective (non-chaotic) bath 
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that is composed of infinitely many oscillators. The proper model that corresponds 
(classically) to the Langevin equation with (|^) is the DLD model that has been 
introduced in [Q. DLD are the initials of "Diffusion Localization and Dissipation". 
These three effects comes out naturally from the quantum-mechanical solution of the 
DLD model. 

Quantum mechanically, it is convenient to use Wigner function p{R, P) in order 
to represent the reduced probability density matrix. We can use the Feynman- Vernon 
(FV) formalism Q| in order to find an expression for the propagator JC{R,P\Ro,Po). 
This propagator is uniquely determined once the friction irf), the temperature of 
the bath (T) and the spatial correlations w{r) are specified. The dependence 
of the propagator on the bath-temperature is via the appearance of a quantum- 
mechanical version of the noise kernel 4>{t). The latter is related to the friction 
coefficient via a universal fluctuation-dissipation relation. At high temperatures the 
quantum-mechanical 0(t) coincides with its classical expression. At the limit of zero 
temperature 4>{t) does not vanish, rather it develops a large negative tail. In the latter 
case, (/'(w), the power spectrum of the noise, reflects the zero point fluctuations (ZPF) 
of the environmental modes. 

The Zwanzig-Caldeira-Leggett (ZCL) model j|] constitutes a special formal limit 
of the DLD model. It is obtained by taking the limit ^ ^ oo. Both the ZCL model 
and the DLD model will be generalized in this paper to allow the analysis of Brownian 
motion in d-dimensions. A modified version of the the DLD model incorporates the 
effect of long range (power-law) spatial correlations. On the other hand, we are not 
discussing the DLD model in its full generality (as in [^), but rather restricting 
ourselves to the particular circumstances that correspond in the classical limit to 
generic Brownian motion. Thus, we are considering the ohmic DLD model, and further 
assume that (||) is satisfied. 

In case of the ZCL model, the quantal propagator IC{R, P\Ro, Pq) is a Gaussian 
stochastic kernel. Consequently the dynamics can be obtained by solving the Langevin 
equation (^ with (^). At low-temperatures the quantum mechanical version of 4>{t) 
should be used. At high temperatures (/>(t) becomes classical-like and consequently 
the ZCL propagator coincides with its classical limit. These observations do not hold 
in case of the DLD model Furthermore, the distinction between the quantal DLD 
propagator and its classical limit persists even in the limit of high temperatures. 

It is important to define what is the meaning of 'high temperatures'. As in the 
classical relatively simple description of the dynamics is obtained if it possesses 

a Markovian property. With a Markovian property it is possible to obtain the long- 
time evolution by composing short-time evolution steps. It is also possible then to 
write down a corresponding Master equation for the Wigner function. The FV path- 
integral expression for the propagator contains the quantum-mechanical version of 
(j){T) rather than the classical one. In order to have a Markovian property one should 
argue that it is possible to use the WNA, meaning to replace (/)(t) by the effective 
classical-like delta-function. Still, in case of the DLD model, the result of the path 
integration is not classical-like: The high-temperature Markovian limit is not the same 
as the classical limit. 

It turns out that the quantum-mechanical condition for using the WNA, thus 
having a Markovian property, is more restrictive than (^ . For ballistic-like motion we 
shall see that the actual condition is, in most circumstances. 

High Temperatures <^ 7 ^ ^ (^) 
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A recent derivation of the high temperature ohmic DLD model which is based on 
a synthetic RMT Hamiltonian has been reported in The existence of such a 
derivation is most encouraging since it further support the idea of having a universal 
description of quantal Brownian motion, at least in the Markovian limit. It has 
been speculated by the authors of the latter reference that a future extension of 
their formalism will lead to an agreement with the general result of The high 
temperature ohmic DLD model has been discussed in Q and has been further analyzed 
in . The physical picture of high-temperature dynamics will be further illuminated 
in the present paper. 

Loss of coherence, or dephasing, is a central issue in the analysis of quantal 
Brownian motion. Wigner's picture of the dynamics leads to the distinction between 
two different mechanisms for loss of coherence. In case of ZCL model one should 
consider the ^spreading mechanism\ This classical-like non-perturbative mechanism 
is very effective in smearing away the interference pattern. In case of the DLD 
model, coherence is much better maintained, and one should consider the perturbative 
^scattering mechanism^ for dephasing. 

At the limit of zero temperature, the spreading mechanism is still effective in 
suppressing interference. One can use the Langevin formalism in order to analyze the 
smearing of the interference pattern. It is important, however, to take into account 
the negative temporal-correlations of the effective noise On the other hand, the 
analysis of low-temperature dephasing in case of the 'scattering mechanism' is a quite 
subtle issue, that constitutes a main concern of the present paper. The Langevin 
formalism is no-longer applicable, and the lack of a Markovian property enforce a 
semiclassical approach. The semiclassical approach has a further advantage: it is 
possible to go beyond the analysis of a simple ballistic-like Brownian motion and to 
analyze other types of transport. We shall distinguish between ballistic, diffusive and 
chaotic motions through cavities. We shall derive various results and contributions to 
the dephasing rate in the various temperature regimes and depending on the physical 
circumstances. Some of these results coincide with similar computations that are 
related to electrons in metal 0. 

An important question is whether the scattering mechanism is still effective in 
suppressing interference at the limit of zero temperature [|[ ^. It turns out that 
our semiclassical approach, in spite of its other advantages, has a limited range of 
validity. It can be trusted if the kinetic energy of the particle is sufficiently large. For 
ballistic-like motion large energy means 

V E 

Large Energy ^ - ^ — (9) 
I n 

in analogy with (^. The latter condition can be cast into the more suggestive form 
A_B <C ^, where As is the De-Broglie wavelength of the particle. The large-energy 
condition is obviously satisfied in case of the ZCL model {£ oo). The large 
energy condition may not be satisfied in case of the DLD model, and consequently the 
semiclassical result should be modified. In particular, in case of a low-temperature 
thermal motion, the contribution of the zero point fluctuations (ZPF) to the dephasing 
rate should be excluded. 
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2. The ohmic DLD Model 

We consider a bath that consists of infinitely many osciUators whose canonical 
coordinates are {Qa,Pa)- The bath Hamiltonian is 

In case of the DLD model the interaction of the particle with the oscillators is described 

by 

Ti-int ^ -y^CaQaM(x-Xa) (ff) 

a 

where Xq, is the location of the a oscillator, m(x— x^) describes the interaction between 
the particle and the a oscillator, and Cq are coupling constants. It is assumed that 
the function u(r) depends only on |r|. The range of the interaction will be denoted by 
i. The oscillators are distributed uniformly all over space. Locally, the distribution of 
their frequencies is ohmic. Namely, 

— ■ S{uj — UJa) S{x — Xa) = fjiO iov iO < \ / Tc . (12) 

It is useful to define the spatial auto-correlation function 

/oo 
u(r-x')u(x')dx' (13) 

Without loss of generality u{y) is normalized such that w"(0) — —1. For example, we 
may consider a Gaussian u(r) for which 

w{v) = £2exp(^-l(£)'^ (14) 

The d-dimensional FT of w{r) will be denoted by w(k). The mode-density (after 
angular integration) is g{k) — {Cd/{'2Tr)'^)k'^^^w{k). In general, we shall assume that 

g{k) ^Cf+'^k"-'^ forfc<l/^ (15) 

where £ characterize the spatial scale of the correlations, and C is a dimensionless 
constant. In case of the short range Gaussian- type correlations (p^), the parameter a 
equals simply to the dimensionality d. For long range power-law interaction it may 
be less than d, possibly negative. In a moment we shall argue that in order to have a 
well defined model we must have [^"(O)! < oo. Therefore only — 2 < cr is meaningful. 
The regime —2 < cr < is well defined but it requires special treatment since w{0) 
diverges. 

The formal limit £ ^ oo corresponds to the ZGL model. The ZCL model describes 
interaction of a particle with environmental modes whose wavelength is much larger 
compared with the distances that are explored by the particle. Consequently, the 
interaction with the a field-mode is approximated by a linear function. The interaction 
term in the Hamiltonian that defines the ZCL model is of the following form (here 
generalized to d-dimensions) : 

'Hint = Vd^CaQaria-X (16) 

a 

The unit vectors fia are assumed to be distributed isotropically. The distribution of 
the oscillators with respect to their frequencies is assumed to be ohmic, as in (^2|) with 
the 5(x— Xq) omitted. 
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The classical analysis of the d-dimensional DLD model constitutes a trivial 
generalization of the one-dimensional DLD model that has been considered in 
Equation (3.10) in the latter reference should be replaced by the following expression 
for the friction: 

/oo 
Q;(r) Vw(vr) dr (17) 
-oo 

where v is the velocity of the particle. For ohmic bath a{T) can be replaced by —r]6'{T). 
One obtains then J- friction = —rj'v. Obviously, in order to obtain a finite result it is 
crucial to have |w"(0)| < oo. Recall that by convention we use the normalization 
w"(0) = —1. Thus, the generalized version dl^ ) of the DLD model is well-defined 
provided — 2 < cr. The reduced motion of the particle will be described by the Langevin 
equation (|l|), where the noise satisfies (||), with 0(r) = 2rikBT5(T). For the realization 
of a classical trajectory, the global functional-form of w{v) is of no importance. Only 
the force correlations (^) are important. The latter correlations are well defined as 
long as |i(;"(0)| < oo, and by our convention they are equal 4>{t). The parameter t 
is insignificant in the classical analysis of the ohmic DLD model, and therefore the 
classical description of Brownian motion is the same as in the case of the ZCL model. 

Quantum mechanically, the reduced dynamics of the system may be described 
by the propagator /C(i?, P|i?o, Pq) of the probability density matrix. For sake of 
comparison with the classical limit one uses Wigner function p(i?, P) in order to 
represent the latter. Using Feynman- Vernon (FV) formalism an exact path- 
integral expression for this propagator is obtained. The FV expression is a double 
sum J J Vx'Vx" over the path variables x'{t) and x"{t). It is convenient to use new 
path variables R = {x'+x")/2 and r — {x"—x'), and to transform the J J VRVr 
integral into the form 

nR,P 

/C(i?,P|i?o,Po) = / T^R IC[R] , (18) 

where IC[R\ is a real functional, which is defined by the expression: 

IC[R] ^ ( Vr e*s-(S/vee+Si.) q-^s« ^ ^^^-^ 



The Vr integration is unrestricted at the endpoints, and the free action functional is 
5/ree[R-, r] = —TO dr Rr. The action functional S'f[R., r] corresponds to the friction, 
and the functional S'Ar[R, r] corresponds to the noise. The explicit expression for the 
friction action functional, assuming an ohmic bath, is 

S'F[R,r] ^ ri I dTVt(;(r)-R (20) 
Jo 

The general expression for the noise functional is 



dtidt2 (t>{t2—ti) X 

[w(x^'-x'/) + w(x^-x;) - 2w(x^'-x;)] (21) 

where is a short notation for x(ti). The power spectrum of the noise (/'(w) is the 
FT of the noise kernel 4>{t). For ohmic bath 

(j){uj) = ^^°^^^('^~f) ^'^^ ^ < l/'^c (22) 
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The ZCL version for Sp and Sj^ is obtained by taking the hmit £ oo which is 
equivalent to the formal substitution w{r) — — r^/2. 

The power spectrum (22) of the quantum mechanical noise is the same as that 
of a classical white noise in the restricted frequency regime lu < ksT/h, where 
^(w) = 2r]kBT. If the temperature is high enough, such that 1/tc < ksT /%, then the 
power spectrum is essentially classical. At lower temperatures the power spectrum 
contains a non-trivial frequency zone kBT/h < uj < 1/tc-, where one may use the 
approximation 0(t) = hrjuj. This component of the power spectrum does not vanish in 
the limit of zero temperature. It corresponds to the "zero point fluctuations" (ZPF) of 
the environmental modes. The temperature becomes effectively zero if ksT/h < 1/t. 
At intermediate temperatures, namely 1/t < ksT/h < 1/tc, it is convenient to write 
the power spectrum as a sum (f>{uj) = 0o(w) + 4't{'^), where 4>a{u}) corresponds to the 
ZPF and (/'t(w) is the excess thermal noise: 

<i>T{Lo) = 2?7 -f^^ for CO < 1/t, (23) 

The propagator (|l8|) possess a Markovian property if it is legitimate to make the 
WNA. The actual condition for having effectively white noise will be discussed 
in a later section. At low temperatures it is essential to consider the non-trivial 
nature the quantal noise in order to obtain the proper quantum-mechanical dynamical 
behavior |10f|. 



3. Propagation in High Temperature Environment 

In the absence of noise and friction, the free-motion propagator of Wigner function is 
the same as in the classical limit. Namely, 

/C(i?,P|i?o,Po) = /CS = 2^,5(P-Fo) '5((i?-i?o) - £<) (24) 
For simplicity, but without loss of generality, we shall consider whenever possible one- 
dimensional motion. In case of the ZCL model, both Sp and Sn are quadratic forms, 
and therefore the propagator is a Gaussian kernel. It follows that the dynamics of 
Wigner function is purely stochastic, and it can be described by Langevin equation 
with (||) and (|2|). In this sense, there are no genuine quantum mechanical effects 
that are associated with Brownian motion, as long as the ZCL model is used for its 
description. The situation is quite different in case of the DLD model. Here the 
propagator is non-Gaussian, and, in general, it may have non-stochastic features. 

In the present section we shall restrict ourselves to the high temperature regime, 
where the white noise approximation (WNA), namely (/)(t) — v5{t) can be applied. 
Whenever the WNA applies, the ZCL propagator becomes identical with the classical 
propagator for damped motion, namely 

/C(i?,P|ii'o,Po) = /cilU = Gaussian (25) 

In contrast, in case of the DLD model, the high temperature limit of the propagator 
does not coincide with its classical limit. An explicit calculation, as in gives the 
following expression: 

JC ^ Wfi/,*/cilUd + e ^\l-Wn/,^)JCn (26) 

This expression is valid also for the generalized DLD model provided < ct. The 
precise value of the parameter £ is defined here via w(0) = P. Above, W{R—R', P—P') 
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Figure 1. Upper plot: Phase space illustration for the structure of the DLD 
propagator. Lower plot: The projected phase-space density. 



is a smooth Gaussian-like kernel that has unit-normalization. Its spread in phase 
space is characterized by the momentum scale h/£, and by an associated spatial scale. 
The symbol * stands for convolution. Thus, the classical propagator is smeared on a 
phase-space scale that correspond to Ap — h/£ and there is an additional un-scattered 
component that decay exponentially and eventually disappears. The structure of the 
propagator is illustrated in Fig.|l|. The significance of this structure will be discussed 
in the next section. 

To write down an explicit expression for the propagator is not very useful. 
Rather, it is more illuminating to follow the standard procedure |^ and to write 
an equivalent Master equation. This is possible since at high temperatures the path- 
integral expression possess a Markovian property. Namely, since both Sp and Sn are 
local functionals of the paths, it is possible to regard the finite-time propagation as 
the convolution of (infinitesimal) short-time kernels. The derivation of the Master 
equation is explained in Appendix B. The final result is 



dp 

at 



Or -P 
m 



Gf*P 



The friction kernel is defined as follows 

and the noise kernel is 
Gn 



(27) 



(28) 
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h/£ 



G 



N 



P-P' 

h/£ 



5{P-P'] 



(29) 



Both Gf and Gn are smooth Gaussian-hke scahng functions, properly normaUzed to 
unity. 

If Wigncr function does not possess fine details on the momentum scale then 
the convolution with Gp can be replaced by multiplication with 1, and the convolution 
with Gjv can be replaced by jdP^ . These replacements are formally legitimate both 
in the classical limit ?i — > 0, and in the ZCL limit i oo. One obtains then the classical 
Fokker Planck equation 



dp 

at 



- Or -P 
m 



dp ^P 



' Qp2 



(30) 



The same observation can be done by inspection of ( p^ . Namely, if the propagator acts 
on a smooth Wigner- function, (no features on momentum scale ti/tj^ then the second 
term becomes vanishingly small, while the first term becomes effectively classical. 



4. Dephasing Within the Wigner Picture 

Wigner function may have some modulation on a fine scale due to an interference 
effect. The standard text-book example of a two slit experiment will be analyzed 
below, where the interference pattern has the momentum scale h/A. See Fig.|[ d is 
the distance between the slits. In this section we shall explain the possible mechanisms 
that lead to the disappearance of such interference pattern. In view of the semiclassical 
point of view of the next section, we use the term "dephasing" in order to refer to 
this disappearance effect. In the present section, as in the previous one, we are still 
limiting ourselves to the high temperature regime. 

In case of the ZCL model, the propagator is the same as the classical one, and 
therefore we may adopt a simple Langevin picture in order to analyze the dephasing 
process. Alternatively, we may regard the dephasing process as arising from Gaussian 
smearing of the interference pattern by the propagator. In case of the DLD model we 
should distinguish between two possible mechanisms for dephasing: 

• Scattering (Perturbative) Mechanism. 

• Spreading (Non-Perturbative) Mechanism. 

Actually, it is better to regard them as mechanisms to maintain coherence. The first 
mechanism to maintain coherence is simply not to be scattered by the environment. 
The second mechanism to maintain coherence is not to be smeared by the propagator. 
The first mechanism is absent in case of the ZCL model. We shall see shortly that 
the first mechanism is much more effective in maintaining coherence. The notions 
perturbative and non-perturbative are used in order to suggest a connection with an 
earlier work |^. The smearing effect can be regarded as arising from accumulation 
of many small- momentum scattering events. A conventional perturbative approach 
cannot be applied in order to analyze such dephasing process, even in the limit of very 
faint noise. 

Having in mind the specific example of a two-slit experiment, we should 
distinguish between two cases. If the environment is characterized by a spatial 
correlation scale £ that is much larger than d, then we can 'forget' about the DLD 
model and just use the ZCL model. The relevant mechanism for destruction of 
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Figure 2. Upper plot: The geometry of a two-slit experiment. The propagation 
of the wavepacket is in the y direction, and the interference pattern is resolved on 
the screen. Lower plot: Phase-space picture of the dynamics. Wigner function of 
the emerging wavepacket is projected onto the {x,px) plane. 

coherence in this case is the spreading mechanism. See detailed analysis in the next 
paragraph. If the environment characterized by a spatial correlation scale £ <ti d then 
a totally different picture emerges. Now, Wigner function contains a modulation on a 
momentum scale much finer than h/£. This modulation is not affected by the friction, 
but its intensity decays exponentially in time. This is based on inspection of either 
the propagator (p^, or the equivalent Master equation (p7|). In the latter case note 
that the convolution with Gn can be replaced by multiplication with —(i/h)'^. The 
decay rate is 

— = -^-4 WNA for < cr 31) 

It is a universal result for the dephasing rate due to the 'scattering mechanism', since 
it does not depend on details of the quantum-mechanical state involved. However, the 
validity of this result is restricted to the high temperature regime, where the white 
noise approximation (WNA) can be applied. 

For completeness we turn back to the detailed analysis of dephasing due to 
spreading. This mechanism leads to non- universal results, since the calculation 
of the dephasing time depends on actual details of the interference pattern. For 
concreteness we consider the simplest case of a two-slit interference experiment. 
The wavepacket that emerges from the two slits is assumed to be a superposition 
V'(^) — S± V'd^l) exp(ifc|x=Fd/2|), where hk is the momentum of the incident particle. 
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and V'(lxl) is a radial envelope function. The Corresponding Wigner function is 

oo 

p(x,p) « (1 + cos(d-p/4)) X pi_,„t(x,p) = P«(x,p) (32) 



n— — oo 

^2/ 



The 1-slit Wigner function is multiplied by cos (id-p^), and therefore it is natural 
to regard the 2-slits Wigner function as composed of partial-wavepackets, each 
characterized by a different transverse momentum. By definition the partial- 
wavepacket p„ equals p for \px — 2Trn/d\ < n/d and equals zero otherwise. Each partial 
wavepacket represents the possibility that the particle, being scattered by the slits, 
had acquired a transverse momentum which is an integer multiplet of Ap^: — 2TTh/d. 
Note that the corresponding angular separation is Ap^/{'hk) = A^/d, as expected, 
the associated spatial separation is Aa; — {Ap/m)-t where t = y/{hk/m) is the 
time up to the screen. It is important to distinguish between the "preparation" zone 
y < d which is excluded from our considerations, and the far-field (Franhoufer) zone 
d^ /Xb <C y where h <C AxApx and consequently the individual partial-wavepackets 
can be resolved. Due to the noise the interference pattern is smeared on a momentum 
scale 6px — V^, and on a corresponding spatial scale Sx = \fvt-t. The interference 
pattern disappears once bpx ~ Ap^ or equivalently bx Ax. This leads to the 
following expression for the dephasing time: 

— = ^^^^d^ WNA for ZCL model, two-slits (33) 

Comparing with (|3l| ) we observe that indeed the smearing mechanism is very effective 
in suppressing the interference pattern. The result is also non-universal since it 
depends on details of the interference pattern. In some other circumstances, where 
Wigner function is characterized by a different type of interference pattern (notably 
the case of "vertical" interference pattern), the dephasing time may be proportional 
to some fractional power of the noise intensity In general, in the full analysis of 
the smearing process, one should take into account the effect of friction. In the above 
example this effect has been neglected. 

The analysis of dephasing in case of the ZCL model is easily extended to the low- 
temperature regime. Langevin formalism is still applicable, provided the appropriate 
4>{t) is being used. The smearing of the interference pattern can be analyzed as in the 
previous paragraph. At the limit of zero temperature one should take into account 
the negative correlations of the noise M. For the ballistic- like motion that has been 
considered above, the smearing scale Sp^ is proportional, at zero temperature, to ln{t) 
rather than to t, leading to an anomalous (non-exponential) dephasing factor M, 0] 



5. Dephasing - The Semiclassical Point of View 

There is a totally different approach to analyze dephasing, that is based on the 
semiclassical point of view. The advantage of the semiclassical approach is the 
ability to pursue the study of dephasing to the low-temperature regime as well as 
to circumstance in which the simple ballistic-like Brownian motion scenario is not 
applicable. Specific examples that will be discussed later are the transmission via 
either chaotic or diffusive cavity. From now on we regard the two-slit experiment as a 
special case of transport problem. See Fig.||. The probability to cross the obstacle, or 
more generally to be transmitted via some cavity, can be written as a (double) sum 
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Figure 3. Various types of transport problems: (a) Ballistic transport as in the 
two-slit experiment; (b) Transport via a chaotic cavity; (c) Transport via diffusive 
cavity, as in weak localization experiments. One should consider also the case of 
ergodic motion via a diffusive cavity (not plotted). 



over classical trajectories. Namely, 

^ A^Al exp ( ,^^WZ^ + _ SM^\ ^^^^ 

ah ^ ' 

where S\yia\ is the classical action for the classical trajectory Xa, and Aa is the 
corresponding classical amplitude. Each pair a&5 of classical trajectories constitutes a 
stationary-phase point of the exact path-integral expression. However, the notation is 
somewhat misleading since once the friction functional is switched on, a and h no longer 
can be considered independent indexes. In particular, strictly speaking, the diagonal 
terms are no longer truly diagonal. However, here comes a very important observation. 
In case of the DLD model, assuming that Xa and x;, are well separated with respect 
to the microscopic scale one has ^'^[xa, Xb] = 0. Thus we conclude that friction has 
no effect on the interference part of the sum (Q) , in consistency with our discussion 
of the scattering-mechanism in the previous section. On the other hand, friction will 
have an effect on the diagonal terms of the sum, as required by the correspondence 
principle. It should be emphasized that in case of the ZCL model, friction may have 
a non-negligible effect also on the interference terms, again in consistency with our 
discussion of the smearing mechanism. 

The suppression of an off-diagonal term in (|3^) is given by the dephasing factor 
exp(— iS'Ar[x£i, Xb]/?i^). Using the definitions of Appendix D it is possible to transform 
expression ( ^ ) into an integral over the Fourier components of the motion involved. 
Consequently, in most cases of interest, the dephasing factor can be cast into the form 
exp(— i/r^), where 

-1 -1 />oo />00 7 

- = -2/ / — 0HP(fc,c.) (35) 

Tip n Jo Jo 

The domain of integration is illustrated in Figj4[ The power spectrum P{k,uj) of the 
motion is calculated in Appendix D for ballistic, chaotic and diffusive motion. Various 
results that can be derived by using the above expression will be presented in later 
sections. 

One should be very careful in the physical interpretation of exp(— i/r^). If we 
have two unrelated trajectories aSzb, and we have also static disorder, then we will 
have "statistical" dephasing that reflects the effect of averaging over realizations of 
the disorder. Consider for example a two-slit experiment: The interference pattern 
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Figure 4. The (fc,aj) plane. Left plot; The shaded regions indicate those 
environmental modes that arc cfTcctivc in the dephasing process. The darker 
region indicates a possible excess contribution due to ZPF. The curve Dk^ 
illustrates the frequency-span of P{k,uj). Right plot: The power-spectrum of 
an ergodic motion via a diffusive cavity is concentrated under the plotted curves. 
See further details in Appendix D. 

on the screen will be distorted for any particular realization of a static disorder, 
and will be washed away completely upon averaging over different realizations. This 
mechanism is non-effective if a&6 are related by time reversal, for which static disorder 
gives S'Ar[Xfj,Xb] = 0. Still some suppression is expected also in the latter case, 
via the classical- like amplitudes |^a|^- In case of dynamical environment the role 
of elastic scattering is taken by inelastic scattering events. Therefore we have genuine 
dephasing, irrespective of whether a&6 are related by time reversal. In some typical 
circumstances, the dephasing factor can be re-interpreted as giving the probability 
to "leave a trace" in the environment (see Appendix C) . We take now the limit of 
zero temperature. Assuming for simplicity that a&6 are related by time reversal, it 
follows from the definition of the influence functional that exp(— 5'Ar[X(j, Xb]/?i^) can 
be interpreted as the probability to excite at least one of the oscillators along the way. 

Up to now we have encountered, depending on the physical circumstances, two 
possible interpretations for the dephasing factor. The "statistical" interpretations 
holds in case of static disorder, while the " leaving-trace" interpretation holds in typical 
cases of a dynamical environment. Let us keep T = and take the limit E ^ 0, where 
E is the available energy of the particle. In order to make the following argumentation 
more illuminating, one may assume that the particle has a small but finite kinetic 
energy E, and that all the oscillators have relatively large frequencies, such that 
E <C uja- Consequently, all the scattering events are elastic. On the other hand 
has a non- vanishing value, which implies that the off-diagonals terms are being 
suppressed. Obviously, the "leaving-trace" interpretation of the dephasing factor no 
longer holds. One wonders whether exp(— S'at/?!^) acquires, under such circumstances, 
a somewhat different physical meaning, as in the case of static disorder. Maybe it 
gives now the probability to be scattered elastically (rather than inelastically) along 
the way, thus leading to dephasing of the "statistical" type. A perturbative treatment 
of the scattering process leads to the conclusion that the probability to be elastically 
scattered by a zero temperature oscillator is proportional to . This is because the 
elastic scattering off an oscillator involves a second-order process of virtual emission 
followed by absorption of a quanta tiUa- At the same time S'Ar[xa,Xh] is proportional 
to to leading order. Therefore, exp(— S'jv/ft^) cannot have the desired physical 
significance. We therefore must conclude that something goes wrong with the present 
semiclassical approach once low energy particle is concerned. We shall come back to 
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this point later on. 



6. Dephasing at High Temperatures 

As a first step in the application of the semiclassical approach, it is interesting to 
recover results (|3l| ) and (^ ) for the dephasing rate. These results hold at high 
temperatures, such that the WNA is applicable. The noise functional in such 
circumstances takes the following simple form 

SnM = iT^ksT f [w{0)-w{r{t'))] dt' (36) 



^0 

where r = x^— Xf,. If the two trajectories are well separated with respect to the 
microscopic scale i, then one indeed recovers the universal result (^). Recall that 
^(O) = In the other extreme limit £ oo, that corresponds to the ZCL model, 
one obtains 

SnM = 2rikBT [ v{t'f dt' = 2r)kBTrl-t (37) 
Jo 

In case of the two-slit experiment r ~ d, and one recovers (|33|). Obviously, it is not 
a universal result. For example, if the separation between the two trajectories grows 
linearly, then the dephasing time will be proportional to 1/T^/^ rather than to 1/T. 

The generalized DLD model with —2 < cr < 0, requires special care since w{0) 
diverges. We are still considering high temperatures, such that the WNA is applicable. 
Using the notations of Appendix D expression (Bq) can be cast into the form 



SnM = i^t / w(r)P(r)dr ut / g{k)dk P{k) (38) 





The latter expression converges for any —2 < cr. In the ZCL limit, where (. is much 
larger than the average transverse distance between the trajectories, one may use the 
small- A; approximation of Appendix D, and then (^) is recovered. In the DLD case, 
where £ is a very small scale, the integration should be split into the domains k < k± 
and k > k±. The wavenumber k± = 1 /r± is associate with the transverse distance r± 
between the two trajectories. Namely, 

{d for ballistic motion t < Tb^i 
^/Dt for diffusive motion r^^i < t < r„g (39) 
L for ergodic motion T„rg < t 

Note that ergodic motion refers to either chaotic or diffusive trajectories with r^,g < t. 
Again we use the notations of Appendix D. For < cr the integration is dominated 
by the upper cutoff k ^ l/£. One recovers then the standard result (|3l|). On the 
other hand, for —2 < cr < 0, the integration is dominated by the k ^ k± modes. The 
wavelength of these modes is like the transverse distance between the trajectories. 
One obtains then 

SnM = 2r]kBT £^ (^-^^ -t WNA for -2 < cr < (40) 

For non-ergodic diffusive motion r± depends on t and therefore the dephasing factor 
is of the type exp{—{t/T^)^^~'^'>^^). In all other cases we have a simple exponential 
dephasing factor. 

Obviously ( ^) is a special case of the general formula (^5|) . Note that ( p5| ) should 
be multiplied by ?i tin order to become a proper expression for Sn- By comparing (BSf) 
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to (|3^) one can reveal the actual condition for the validity of the WNA. Within the 
interval < k < l/£, the power spectrum of the motion occupies a frequency range 
Lv < LOch where the effective cutoff is 

_ J v/£ if A;=l/£ is located within the ballistic regime 
'^^ \D/£^ if A;=1/£ is located within the diffusive regime 

As long as ujci < ksT/h one can use the WNA in estimating the integral in (^sf). 
One observes that under such circumstances ( ^5|) reduces, after multiplication by fi t, 

7. Dephasing at Zero Temperatures 

As the temperature becomes low, such that ksT /fi < lucI, the dephasing rate becomes 
larger than the value which is predicted by the WNA. This is due to the "zero point 
fluctuations" (ZPF) in the frequency zone kBT/h < uj. The temperature becomes 
effectively zero if ksT/h < 1/t. Time longer than t is required in order to resolve 
such low temperatures. One can use then the zero temperature limit of 0(r) for which 
(j)(Lu)=hr]LU. This power spectrum corresponds to the ZPF of the environmental modes. 

For either ballistic or chaotic motion. The integral ( p5| ) is dominated by {k,uj) 
modes that are concentrated below the curve uj — vk. One obtains 
1 CI 

— « — X —Ti£ V for ballistic or chaotic motion (42) 

(l+cr)7r n' ^ ^ 

Similarly, for diffusive motion the integral is dominated by {k,uj) modes that are 
concentrated below the curve u — Dk^ . One obtains 

1 C , / ftvV\ 1 



, In 1+ — X —riD for diffusive motion (43) 
(2+a)7r y \DJ j n' ^ ' 

In both cases most of the contribution comes from modes with large wavenumber, 
namely k^l/£. 

The validity of the present semiclassical approach, which is based on the 
stationary-phase approximation is limited. Common wisdom [ pT| is that the 
applicability of the semiclassical approach is restricted to circumstances such that the 
energy transfer between the particle and the environment is much smaller than the 
particle's available energy. Technically, it is equivalent to the assumption of uncoupled 
wave-equations [Q. The coupled wave-equations for the particle-oscillators system 
can be uncoupled provided certain conditions are satisfied. Then we can treat the 
particle as moving with constant velocity v and solve for the oscillators. It turns out 
that this reduction requires the assumption of small energy transfer. Therefore, one 
should anticipate problems once oscillators with larger than E are involved. In 
the latter case, there is no justification to think of the particle as decoupled from the 
bath, moving with some constant velocity, capable of exciting oscillators along the 
way. Therefore the corresponding factor exp(— S'at/Ti^) loses its physical significance. 
From the above considerations it follows that a reasonable condition for the validity 
of our semiclassical approach is 

hcuci < E definition of large energy (44) 

On physical grounds this is the condition for being able to leave a trace "along the 
way" . The semiclassical significance of this condition can be further illuminated. 
For ballistic-like motion it is equivalent to the condition of small momentum transfer 
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h/i < p. This is precisely the condition for the applicability of the semiclassical 
methods for the scattering of the particle by the oscillators. More generally, one 
can define a quantum mechanical version of P{k,uj) as the FT of the correlator of 
the Heisenberg-picture operator exp(zfcx(i)). The quantal P{k,uj) coincides with the 
classical P{k,uj) only for uj < E/h. Larger classical frequencies are not supported 
by the effectively banded energy spectrum. Our practical conclusion is that the 
contribution of ZPF should be included if non-thermal motion with large energy is 
concerned. It should be excluded in case of thermal motion with E ^ ksT. 



8. Dephasing due to Thermal Noise 

The possible contribution of ZPF to the dephasing has been discussed in the previous 
section. Now we are interested in the the thermal noise contribution (TNC). The TNC 



can be calculated using (35) with 0(w) replaced by 4>t{^)- See ( |23| ) for the definition 
of the latter. As the temperatures is being raised various regimes are encountered: 

Extremely Low Temperatures <;4> 1/t <kBT/h<l/T± 

Low Temperatures <4> 1/t± < ksT /fi < lOci 

High Temperatures ^ LUd < ksT/fi (45) 

The time scale tj^ corresponds to the transverse distance r±_ of (^9[). It equals d/w for 
ballistic motion, and t^^^ for ergodic motion. In case of non-ergodic diffusive motion, 
where t_l = the extremely low temperature regime is absent. At high temperatures 
we can use the WNA as discussed in previous sections. The low and the extremely 
low temperature regimes are discussed in the next two paragraphs. 

When considering extremely low temperatures it is useful to define the critical 
exponent CTc = 1 for ballistic motion, and ac — 2 for diffusive motion, li a < Uc the 
integration in ( |35| ) is dominated by fc ^ fc^ and one obtains the result 

1 _ ^,2T^kBT fey fkBT/h 



For (Tc < f the integration is dominated by fc ~ 1/^ and one obtains the result 

1 _ ^, X (47) 



In the formulas above C is a numerical factor of order unity. Note that in both cases, 
disregarding possible contribution of ZPF, the dephasing rate is proportional to T^. 
This should be contrasted with the high temperature behavior where, disregarding the 
special case of non-ergodic diffusive motion, the dephasing rate is proportional to T. 

At low (but not extremely low) temperatures it is useful to define a temperature 
dependent wavenumber as follows: 

7 _ / {kBT/hv) in the ballistic regime . 
" 1 (kBT/hD)^/^ in the diffusive regime ^ > 



For (T < the integration in (^5|) is dominated by fc ~ /cx, and one can use the WNA 
result (|40|). Thus, for cr < 0, in the absence of ZPF contribution, we can trust the 
WNA at both high and low temperatures, and an actual crossover is expected only 
when extremely low temperatures are involved. For < a < ac the integration is 
dominated by fc ~ fc^, and one obtains the result 
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If (Tc < a the integration is dominated by k ^ l/l and one obtains again (^7|). 
Thus, for < (T < tJc, there is a non-trivial low-temperature regime where the 
TNC to the dephasing rate is proportional to a non-universal power of T . Below 
this intermediate temperature-regime the TNC is proportional to . Above this 
intermediate temperature-regime we can trust the WNA and the dephasing rate is 
proportional to T. 

9. Manifestation of Effective Static Disorder 

We shall consider in this section quantal Brownian motion in 2 < d dimensions, 
that is described by the the DLD model with short range spatial correlations. For 
concreteness let us assume that the fluctuations of the effective stochastic potential 
satisfy: 

(Z^(x",t")Z^(x',t')> - 0(t"-t') -^'exp l^-i (^^)'^ (50) 

If we are interested in the dynamics over a finite time interval i, then obviously all 
the Fourier components in the frequency regime |a;| < 1/t will have the same effect 
as static disorder. Let us denote by U{x,t) the "static" component of the effective 
stochastic potential. The variance of U is determined by the product (j>{uj~Q) ■ {1/t). 
Recall that (/)(cl'=0) = Consequently we have 

{Ui^"X)U{^\t')) - exp ( ( ^^-^ ] ) (51) 




where = vt"^ /t. We would like to find out whether this effective static disorder 
will manifest itself. Similar question has been raised in ||l^. One should not take for 
granted that the effect of low-frequency fluctuations is completely masked by the the 
incoherent effect of the high frequency modes. 

The first obvious step is to calculate the "statistical" dephasing rate due to the 
presence of the (effective) disorder. We can use (^) with the formal substitution 
(/)(t) = thus obtaining 

■ / g{k)dk P{k,uj=Q) (52) 
Jo 



T, 



^ P JO 



For ballistic-like motion one obtains 

V \n\ t 



n - -^W^- = ( ^ ) • - (53) 



In order to have a non- vanishing effect we should have <C i. This condition can be 
cast into the form ^ <C where ^ = {fiv)'^ / {£W'^) is the mean free path. Still another 



way to express this condition is r™"^ <^ i/v, where r™"* is given by (|3l|). On the 
over hand, the actual dephasing time should satisfy £/v < r^, else coherent effects due 
to the scattering by the (effective) disordered potential will not manifest themselves. 
Thus, we come to the conclusion that the following condition should be satisfied in 
order to have manifestation of coherent effects due to scattering by the effective static 
disorder: 

r*~- « - « (54) 
Obviously, this condition can be satisfied only at the low temperature regime. 
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10. Dephasing versus Dissipation, Concluding Remarks 

It is quite striking that the friction coefficient i] is not affected by quantum mechanical 
effects. Having the Boltzmann picture in mind it is not anticipated that aU the 
quantum mechanical scattering events will conspire to give the classical result. Still, 
this conclusion follows from the FV-formalism quite easily. In order to have a clear 
physical picture, let us consider the time evolution of a Gaussian wavepacket using 
either the propagator (|2^) or the equivalent Master equation (|2^). Clearly, the 
dissipation rate dE/dt is the same as in the classical picture, with small transient 
corrections. Thus, also in the quantum-mechanical picture, the damping process is 
characterized by the time constant = (r]/m)~^. 

Consider a ballistic- like motion with non-thermal energy, such that h/£ <C p. 
Consequently the equivalent condition is satisfied, and we can trust semiclassical 
considerations as far as dephasing is concerned. We can cast the universal WNA 
result (|l]) and the ZPF result (H) into the form 

1/Ty f {£/\Tf high-T (universal) , . 

1/t„ ~ 1 (i/XB) low-T 

where At is the thermal wavelength, and As is the De-Broglie wavelength. We see 
that for the above non-thermal motion the dephasing time is always shorter than the 
damping time. The dephasing rate is linear in T at high temperature, and saturates 
at low temperatures. 

For thermal motion the latter statement is no-longer true. At the low temperature 
regime the ZPF contribution should be excluded, and therefore the dephasing rate goes 
to zero, while the damping time remains finite. One should distinguish between various 
physical circumstances. At extremely low temperatures, assuming ergodic motion, the 
dephasing rate is proportional to T^. At higher temperatures we can trust the WNA 
provided ct < 0. Otherwise, there is an intermediate low-temperature regime where 
the dephasing rate is proportional to a non-trivial power of T. 

In this paper we have introduced a systematic derivation of a general formula 
for the dephasing rate, Eq.(|35|), that holds in all physical circumstances, including 
the ZCL limit. Expressions similar to Eq.([3^), that incorporate integration over the 
(fc, to) environmental modes, are encountered frequently in the literature, starting from 
the well known work in JT4[ . Whenever diffusive motion of electrons is concerned 
(see Appendix E), our expressions agree with well know results |7|. However, most 
publications avoid a straightforward application of the FV formalism, and introduce in 
some stage heuristic considerations in order to obtain a convergent result. There is a 
"zoo" of cutoffs that are introduced in performing the (fc, w) integration, some of them 
are questionable. For example, it is customary to take l/r^, itself as a lower cutoff 
for the uj integration. The upper cutoff is sometimes ksT/h, sometimes Fermi energy, 
sometimes the kinetic energy that corresponds to the drift velocity, sometimes the 
inelastic scattering rate, and so on. The role of Fermi statistics in the determination 
of the various cutoffs is usually left unclear. There are similar ambiguities in the 
determination of the the proper cutoffs for the k integration. Our derivation has 
led to a proper definition for the power-spectrum of the motion, and has made it 
unnecessary to introduce ad-hoc cutoffs into the calculations. Furthermore, in the 
analysis of dephasing at the limit of zero temperature, we were successful in reducing 
our considerations to the level of "one particle physics", thus avoiding a complicated 
discussion of the role played by Pauli exclusion principle. 
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Appendix A. Useful Identity 



We first cite the well known identity 
1 



: exp 



1 



7:Stdl 
2 



0{St^) 



5{x - xo) 



V2nSt V 25i 

Both sides of this identity should be interpreted as kernels of operators. When applied 
to wave-functions, the left hand side corresponds to the free-motion propagator, and 
the second term in the right hand side corresponds to the free-motion Hamiltonian. 
We shall explain now how to derive the following related identity: 



1 



■ exp 



-{R- Rn){r-rn)-iijiR- Ro) 
dt 



2TrSt 

[1 + iStdndr + ijStdr + 0(6t^)] S{R - Ro)S{r - ro) (A.l) 

Both sides of the latter identity should be interpreted as kernels of operators that 
operate on phase-space functions. We start the derivation by writing an equality that 
follows from the first identity via simple replacements: 



, : exp ± ^ , 

Upon multiplication it follows that 

1 ■ * / 2 

;{X^-X 



1 



-Stdi + OiSt'') 



S{x± Tiv^t) 



i 

1 + - 



exp ± 



6tdl 



2St 

S{x+) 



.1 . 

i—ri[x^ 



1 



r]5td'{ 



Simplification of the right hand side gives 

i + Lst{dl-d'_) + ^vSt{d, 



5{x_) + -T^5t5'{x+) 



6{x+)5{x^) 



The replacements x+ —t (x" — Xq) and a;_ (x' — x'q), followed by the transformation 
to the variables R = {x" + x')/2 and r = (x" — x'), gives the desired result. Note that 
the derivation holds also in the case where t] is replaced by some function of r. 



Appendix B. Derivation of the Master Equation 

Wigner function p{R, P) is the Fourier transform of the reduced probability density 
matrix p{R^ r) in the variable r ^ P. The path integral expression for the kernel 
IC{R,r\RQ,ro) is the type J jDRDr with an obvious endpoint conditions, and its 
integrand is essentially the same as in (|l^). See ||^ for more details. The infinitesimal- 
time kernel equals simply to the integrand of the path integral expression, namely 



ICst - ciip(--{R-Ro){r-ro)+i^w'{r)iR^Ro)~dt-^{w{0)-w{r)) 
^ dt n m h 



1 + St( i-dndr - -w'Mdr - ^(w(0) - w(r))] + 0(6t^) 
\ m m h J 



S{R~ Ra)S{r-ro) 
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where the second hne is obtained by employing the identity of Appendix A. Thus, 
to leading order, the kernel ICst has the same effect as operating with a differential 
operator of the type 1 + StC. Consequently, we find that p{R, r) satisfies a Master 
equation of the form dp/dt = Cp. This Master equation is easily transformed into 
the Wigner representation, where any derivative dr is replaced by multiplication with 
iP/Ti. Similarly, any multiplication by a (real symmetric) function G(r) is transformed 
into a convolution with a (real symmetric) kernel G{P—P'). For the convolution 
/ G{P-P')p{P')dP' we use the notation G ★ p. A kernel G{P-P') and a function 
G{r) are related by a Fourier transform with the convention 

f / P-P' \ 

In the above expression we have modified the standard FT convention, by including 
the factor l/(27r?i). This has been done in order to have properly normalized kernels, 
with the measure dP rather than dP/{2TTh). Note also that the FT ofw'{r) equals to 
the iUdp derivative of the kernel Gp that corresponds to the real symmetric function 

w'{r)/r. 



Appendix C. Dephasing and Inelastic Scattering 



Consider the interference contribution of two trajectories a&6 that are related by time 
reversal. Assuming short range interaction with the environmental modes we can prove 
that dephasing is related to the probability for leaving a trace in the environment. This 
statement is true at any temperature. "Short range interaction" means that < a 
and that f is a small scale. " Leaving a trace" means that at least one of the oscillators 
has changed its quantum-mechanical state. It follows that under such circumstances 
dephasing-rate is equal to the inelastic scattering rate. 

Recall the definition of the influence functional. For simplicity we consider first a 
zero-temperature bath, meaning that all the bath-oscillators are initially in the ground 
state. These oscillators are driven by the motion of the particle. The excited states 
of the bath will be denoted by The evolution operator of a driven oscillator 

will be denoted by Ua[x]. 



F\Xa,Xb] 




{na} 



X6 



{0} M 



{0} 



If aSzb are related by time reversal, then one can consider only those oscillators that 
are located "along the loop". Each oscillator along the loop is characterized by its 
natural frequency LUa and by the time ta at which Xa{t) ^ Xa- Consequently we have 
the following expression for the influence functional: 



F[xa,Xfc] = e 
where the excitation probabihty is F({nQ}|{0}) = H K^-al^al^allO)!^ 



(C.l) 



We would 



like to argue that F[xa,Xf,] — P({0}|{0}). Indeed, the summation in (C.l) contains 
one-oscillator excitations for which — 1, two-oscillator excitations for which 

^ Uq, = 2 and so on. Let us consider all the one-oscillator excitations that involves cJq, 
in the range [uj,u!-\-6uj]. Each of these excitations contributes the same P({nct}|{0}), 
but with a different phase factor exjp(i2ujata). By construction of the DLD model 
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(Eq.(12)), the summation over the phase factors will lead to a zero contribution. The 
meaning of "zero contribution" is as follows: either we average over realizations of Xa, 
or else we recall that the influence functional appears inside a path-integral expression. 
The summation over paths will have the same effect like an averaging procedure, and 
therefore we will have indeed a " zero contribution" in the mathematically-generalized 
sense. Similar argumentation applies to other subsets of excitations. If the bath is 
initially in thermal equilibrium, rather than in zero temperature, than the relation 
F[xa,Xh] — -P({0}|{0}) can be generalized into 

F[xa,Xft] = exp(— S'Ar[xQ, Xfc]/?i^) = P(leaving no trace) . (C.2) 

In order to prove the latter statement one should use the same procedure as above with 
|{0}) replaced by some arbitrary initial preparation Then the result should 

be thermally averaged. 



Appendix D. Statistical Characterization of the Trajectories 

Having a pair a&6 of trajectories, we can define the function 

Pabir, t) - (^(r - (xjt'+r) - Xfc(t')))) (D.l) 

where the average is over t' within the time interval [0, t] which is considered. We 
shall use the notations Paa ~ P\\ and Pab = P± for a ^ b. It is also useful to define 
the functions 

P(r,T) = P„(r,r)-P^(r,r) 

P(r) = P(r,0). (D.2) 

The Fourier transform of P(r) will be denoted by P(k), and the double Fourier 
transform of P(r, r) will be denoted by P(k, oj). The above functions all appear within 
integrals where "isotropic" integration over r or k is being performed. Therefore 
it is convenient to average all these functions over all orientations, and so to have 
functions that depend on either r = |r| or fc = |k| respectively. A useful notation 
is Cos(r) = (cos(J7-r)), where the average is over the orientation of the unit vector 
fl. The function Cos(r) depends only on |r|. It equals cos(r) in 1-D, regular Bessel 
function J^o(r) in 2-D, and sinc(r) in 3-D. 

We shall distinguish now between ballistic trajectories as in Fig.^(a), chaotic 
trajectories as in Fig.^(b), and diffusive trajectories as in Fig.||(c). We can also add 
to this list the case of diffusive trajectories that cover ergodically the whole available 
space. In case of a ballistic trajectory P||(r,r) — S{r—r), leading to i^, (k, r) — Cos{kf), 
where f — vt. Similar expression hold for P^, where f « (d^ + (vr)^)^^/^, and d is 
the transverse distance between the slits. For a chaotic trajectory Pj (r, r) starts as in 
the ballistic case, while Pi(r, t) is uniform in r. The ergodic time is given essentially 
by the ballistic time, namely T„g = Ttai = L/v, where L is the linear dimension of 
the cavity. For re,.g < r both P|| (r,T) and P_i(r, r) become uniform in r. For diffusive 
motion P||(r,T) is as for the ballistic case as long as r < Tb^i where Tbai — D/v'^. On 
larger times it becomes a Gaussian, and consequently P||(k, t) = exp(— Ufc^r), where 
D is the diffusion coefficient. The latter expression as well as the approximation 
Pi (r, t) k, (r, t) hold for non-ergodic diffusive motion. After the ergodic time 
T„g = I? jD^ both P||(r, r) and Pi(r, t) become uniform in r. 

In the following paragraphs we discuss the various fc regimes of the function 
P{k,uj). We shall distinguish between the small-fc regime and the large-fc regime. In 
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case of diffusive trajectories the large-fc regime wih be further divided into a baUistic 
regime {v/D < k) and a diffusive regime (fc < v/D), where v/D is the inverse of the 
mean free path. We turn first to define the notions of small and large k. The transverse 
distance between two trajectories has a distribution P±{r), and we can define a typical 
value r± = ((r^))^/^. See Eg. (|39|) . The associated wavenumber is k± = l/r±. Large 
k means from now on k± <ti k. For large k we may use the approximation 

P{k,uj) « P||(fc,CL') for large fc (meaning fc^ <C fc) (D.3) 

For either ballistic, diffusive or chaotic trajectories we have 

P{k,uj) K (■^\ in the ballistic regime (D.4) 



vk \vk 

where B is the FT of the Cos function. This rectangular-like scaling-function has a 
unit width and a unit normalization. In case of diffusive motion we also have 

2Dk^ 

P(k,uj) w , „, in the diffusive regime (D.5) 

Diffusive regime means here any large k that satisfies k < v/D. The validity of the 
latter approximation is further restricted by the condition uj < v^ / D. The collision 
frequency oj ^ / D should be used as a cutoff to the slow power-law decay. In 
our calculations we shall assume that l/£ ^ v/D, meaning that the motion is diffusive 
also on the spatial scale 

For small fc, meaning k < k±, the transverse term cannot be ignored, and we no- 



longer can use the approximation (D.3). We can define a time scale t± that correspond 
to k±. It equals d/u for ballistic motion, t for diffusive motion, and r„g for ergodic 
motion. For a given small k, the function P{k, t) is concentrated within t < t±, since 
i-^i and P± are not identical there. Consequently 

P{k,uj) « T±-P{k) for UJ < 1/t± in the small k regime (D.6) 

The function /^(k) is the FT of -Pj_(r), and therefore we have for P(fc) = I— ^^(k) 
the following small k approximation 



kd^k^ for ballistic motion 



P{k) w ^^i^^ = \ Dt-k'^ for diffusive motion (D.7) 
\L^k^ for ergodic motion 



For large k we have P(k) « 1, as implied by (D.3) and the definition (p.2|) 



Appendix E. Friction Constant For Electrons in Metal 

The effect of electron-electron Coulomb interaction can be analyzed by considering 
the motion of a single electron under the influence of a fluctuating electrostatic 
potential U{x,t) that is created by the other electrons. Thus the electron experience 
an interaction with a fluctuating field that is characterized by {U{x" ,t")U{x' ,t')) . It 
is well known, using fluctuation-dissipation theorem, that for diffusive electrons the 
corresponding fluctuation spectrum is 

w(k)(j){uj) = — f^TT"^) forw<— , \k\<l 
ad k'' \2kBi J Tc t 

The ohmic behavior is cut-off by the Drude collision frequency I/t^. The elastic mean 
free path is £ — vTc, where v is the Fermi velocity. For k < l/£ the power spectrum 
P(k,uj) of the diffusive motion is concentrated below the Drude cutoff frequency. 
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The ballistic regime l/£ < k is oi no interest because its contribution is suppressed 
due to the Drude cutoff. One observes that these ffuctuations corresponds to the 
generalized DLD model (|l5|) with a = d—2. Having observed that the mean free path 
is the physical (effective) cutoff for the spatial fluctuations (this statement is true for 
(T<4), it follows from the convention |z/;"(0)| = 1 that the friction constant (up to 
dimensionless factor of order unity) is given by the expression 

~ iV 1 



Here e is the charge of an electron, and ad is the conductivity defined for d — 1,2,3 
dimensions. The last equality is obtained by using Einstein relation in order to express 
ad in terms of D. The notation stand for the mean level spacing within a cube 
whose volume is i*^. 
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